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ABSTRACT 

We analysed the archival XMM-Newton data of radio loud narrow line Seyfertl 
galaxy PMN J0948+0022 in the energy range 0.3-10.0 keV. The X-ray data 
revealed the spectrum in 0.3-10.0 keV energy band is not a simple power-law 
as described in literature. Instead it consisted of a power-law with soft excess 
below 2.5 keV. The X-ray spectrum was fitted with four different models and 
it was shown the soft excess component of the spectrum in 0.3-2.5 keV energy 
range could be described reasonably well within the framework of the thermal 
Comptonization model as well as relativistically blurred reflection model. The 
power-law component required to fit the spectrum beyond 2.5 keV was found 
to be rather hard compared to the ones observed in other obscured active 
galactic nuclei. It is also shown that the Swift/XRT spectrum from the source 
could not reveal the soft excess component due to poor statistics. The frac- 
tional variability estimated from XMM-Newton data indicates the presence 
of independently varying components in the spectrum above 1 keV. 
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1 INTRODUCTION 



PMN J0948-I-0022 is a radio-loud narrow line Seyfert 1 (RLNLS l) galaxy, with cosmological 



redshift z = 0.5846 and radio loudness parameter R ^ 846 (IFoschini 



7-rays by Large Area Telescope (LAT) on board Fermi satellite (lAbdo et al.l 



2011). detected in 



2009b 



also one of the six such sources detected by Fermi/LAT in the MeV-GeV regio: 



I 



a). It is 



Abdo et al 
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(j2009bl Jal) report ed the multi-w a yeleng th observations of PMN J0948+0022 during March- 



July 2009, while 



Foschini et al. 



(120111 ) observed the source during a 7-ray outburst in July 
2010. In both cases the observed spectral energy distribution (SED) looks very similar to 
the SED of typical high power blazars, known as flat spectrum radio quasars (FSRQs). 
The SEDs were model e d wit h the synchrotron and Comptonization models developed by 
Ghisellini fc Tavecchio ( 2009I ). For both the observations the broadband SED were fitted 
with a synchrotron component for radio-to-IR/optical, a synchrotron self-Compton com- 
ponent for the X-rays, external Comptonization component for MeV-GeV 7-rays, and a 
thermal component from the accretion disk in the UV band. Interestingly, the fitted param- 
eters and estimated power were found to fall in the region between the low power blazars 
(BL Lacs) and high power blaza rs (FSRQs). These derived properties are similar for other 
three Fermi observed RLNLSl (lAbdo et al.ll2009a] ). Their analysis of broadband SED lead 
to similar conclusions. But the most intriguing fact is that the estimated black hole (BH) 
mass and mass accretion rates for RLNLSl, in terms of Eddington rate, are very differ- 
ent compared to those properties of blazars. The BH mass for RLNLSl is approximately 
two orders of magnitude lower than the BH mass in blazars (Foschini 2011). Secondly, the 
RLNLSls accrete at near Eddington limit, while blazars accrete at a much lower rate. 
Here in this work, we gave a closer look at the X-ray spectrum of PMN J0948+0022. The 
reason is the following. In the broadband SED modelling for all 7-ray loud RLNLSl the 



X-ray component is described by synchrotron self-Compton emission process (lAbdo et al 



2009b 



a|). The X-ray spectra of radio quiet narrow line Seyfert 1 galaxies, as observed with 



XMM-Newton , show the presence of soft excess in the 0.3-2.0 keV energy band and a power 



law componen 


t ( 


Bian & Zhao 


2004; 


Boiler et al. 


2003 




i^'abian et al. 


2004 


O'Brien et al. 


2001; 


Page et al. 


200,4 


Dewangan et al. 


2007 


). The origin 



Haba et al. 



2008 



is not properly understood yet. In fact some of radio emitting NLSl which are not catego- 



rized as radio-loud, a . 



504 (iPapadakis et al. 



so sh ow strong soft excess in their X-ray spectrum, e.g, P KS 0558- 



20101). PKS 2004-447, another radio -loud (R = 6358) NLSl ((Foschini 



2OIII ). also show strong soft excess in its X-ray spectrum (IGallo et aL 



2006). But there are 



some doubts about the classification of the source PKS 2004-447 (jOshlack et al. 



20011 ). If 



the X-ray spectrum of 7-ray loud RLNLSl PMN J0948+0022 also show the presence of soft 
excess then it is important to understand their origin as well as their connection to the 7-ray 
emission. Here, we discussed the analysis and modelling of archival data of PMN J0948+0022 
observed with XMM-Newton , Swift/XBT to show that the X-ray spectrum of the source has 



The XMM-Newton view of RLNLSl galaxy PMN J 0948+ 0022 3 

a complex structure and Swift/XBT can not adequately describe it. We also analysed the 
Fermi/LAT data taken simultaneously with XMM-Newton observation to study the 7-ray 
emission from PMN J0948+0022. 

The paper is structured as follows. In section 2, we described the data reduction for XMM- 
Newton, Swift/XBT and Fermi. The spectral modelling is described in section 3, results are 
discussed and concluded in sections 4 and 5, respectively. 



2 DATA REDUCTION 



2.1 XMM-Newton observations 

PMN J 0948+0022 was obs erved by XMM-Newton telescope (for details see 



Striider et al 



( I2OOII ): [Turner et al.l ( 1200 ll )) operated with Medium filter at two epochs during April, 2008 
(Observation ID 0502061001, hereafter Obs-A) and May, 2011 (Observation ID 0673730101, 
hereafter Obs-B) for approximately 24 ks and 93 ks, respectively. The X-ra y data is not 
taken in prime full window mode for EPIC-MOS detectors (MOSl and M0S2; 



Turner et al. 



( I2OOII )) for Observation ID 0502061001 and the source was not ava ilable within the CC D 



frame, therefore, only EPIC-PN (European Photo Imaging Camera; IStriider et al.l (120011 )) 
data were used for further analysis. For Obs-B, the total observation time intervals of MOS 
detectors were split into five short duration segments of 16 ks (Ml), 28 ks (M2), 25 ks (M3), 
4 ks (M4), 10 ks (M5), however, the data have been taken for full observation time span, 
i.e. 93 ks, for PN detector. 



2.1.1 EPIC data 

Data reduction followed standard procedures using XMM-Newton Science Analysis System 
software (SAS version SASll.0.0) with updated calibration files. Event files for the MOS 
(for Obs-B) and the PN (for Obs-A and Obs-B) detectors were generated using the tasks 
Emchain and Epchain, respectively. These tasks allow calibration of the energy and the 
astrometry of the events registered in each CCD chip and to combine them in a single data 
file. Event list files were extracted using the SAS task Evselect. Data from the three 
cameras were individually screened for the time intervals with high background when the 
total count rate (for single events of energy above 10 keV) in the instruments exceed 0.35 
and 1.0 counts s^^ for the MOS and PN detectors respectively. 

Source photons were extracted from circular region of 40" around the source position to 
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generate the light curves and spectra. Background was estimated from a number of empty 
regions close to the X-ray source in the detector. X-ray spectra and light curves of the 
source were generated using SAS task Evselect and background scaling factor have been 
calculated using Backscale task. The photon redistribution as well as the ancillary matrix 
were computed using SAS task Rmfgen and Arfgen. Finally, the FN spectra of Obs-A 
and Obs-B were rebinned to have at least 20 and 80 counts per bin respectively while the 
MOSl and M0S2 spectra of Obs-B were rebinned to have at least 20 counts per spectral 
bin. 



2.2 Swift observations 

X-ray data from X-ray telescope (XRT) on Swift have also been analysed here in photon 
counting mode. These observations have been taken on the same date as for Observation ID 
0673730101 from XMM-Newton, i.e., 28 May 2011. The Swift/XRT data has been reduced 
using the HEAsoft packages ( 6.12 version), including Xselect, Ximage, and Swift data 
analysis tools. The XRT tool Xrtpipeline is used to produce the clean event data. X- 
ray spectra were extracted using Xselect with 20 pixel radius source extraction region 
and background region from clean event data. Response matrices files (RMF) and Ancillary 
Response Files (ARFs) were generated from task Xrtmkarf using the Caldb version 
014 (released date 25 July 2011). Finally, the spectra were grouped using the Grppha task 
with 20 counts per bin. 



2.3 Fermi/LAT observation 

The Large Ar ea Telescope(LAT), o nboard Fermi gamma-ray observatory, is a pair conver- 
sion detector (lAtwood et al.ll2009l ). It is sensitive to the gamma-rays in the energy range 
~20 MeV to >300 GeV with peak effective area of ~8000 cm^ at 1 GeV. LAT has a very 
wide field of view of 2.4 sr covering 20% of the sky at any instant with angular resolution of 
~1 deg. at 1 GeV. 

The LAT data for the source considered for the present analysis cover the period May 15- 
June 15, 2011. The data have been analyzed using the standard Fermi-L AT ScienceTools 
software package (ver. v9r27pl)^ and with the help of the analysis threads and other doc- 



^ http:/ /fermi. gsfc.nasa.gov/ssc/data/analysis/software/ 



..tit. 
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OBS ID 0502061001 fPN) : Power law 




(a) 



MOS1 : Power law model 



mmm 



mm' 



(c) 



Power law model 
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Energy (keV) 

(b) 

MOS2 : Power law model 




(d) 



Figure 1. Soft excess is shown for (a) EPIC-PN spectrum for Obs-A, (b) EPIC-PN, (c) MOSl and (d) MOS2 spectra for 
Obs-B. The power-law model is fitted for the energy range 2.5 keV and above in each case. 



umentation available at Fermi Science Support Center webpageqj. We have selected only 
"Diffuse" class events (i.e. the purest gamma-ray events) in the energy range 100 MeV 
to 300 GeV within 10° region of interest (ROI) centered around the source position (RA 
= 09^8*^57^3 and Dec = +00°22'25".6) using gtselect task. In order to avoid back- 
ground contamination from the bright Earth limb, we have discarded photons arriving from 
zenith angle greater than 105°. We selected only those time intervals during which the 
events/photons qualify for science analysis using the filter expression (DATA_QUAL==1 
kk LAT_C0NFIG==1 kk ABS(ROCK_ANGLE)<52) within gtmktime task. Then the 
livetime for the source was calculated and exposure map of the source was generated. Fi- 
nally, using Unbinned likelihood analysis tool (gtlike), we fitted the source spectrum along 
with other 21 sources present within the ROI of 10° centered around PMN J0948+0022 and 
obtained the flux in the energy band 100 MeV - 300 GeV. For this we created source model 



^ http:/ /fermi. gsfc.nasa.gov/ssc/data/analysis/ 
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Multicolour disc blacl<body and power law model 



Thermal Comptonization and power law model 




(a) 

Smeared wind absorption model 




(c) 




(b) 

Belativistically blurred rellection model 




(d) 



Figure 2. EPIC-PN X-ray spectrum of Obs-B with best fit models : (a) Multicolour blackbody with power-law, (b) Thermal 
Comptonization and power-law, (c) Smeared wind absorption model and (d) Relativistically blurred reflection model, 
distributions of fitting are given in lower panels of each plot. 



file containing model spectra for all these 21 gamma-ray sources visible in the ROI as listed 
in the 2nd Fermi/LAT catalogue and diffuse emission models of galactic and extragalactic 
origins. For PMN J0948+0022, we have taken power-law spectrum model. We obtained the 
value of test statistics (TS; detection significance of the source ~ VTS) as 32.86 for the 
whole energy range 100 MeV-300 GeV. 



3 RESULTS 

3.1 Spectral analysis 

The EPIC-PN spectra of both the observations and MOSl and M0S2 spectra of Obs-B were 
fitted with a power-law dis tribution with the photon index F using the spectral modelling 



tool XsPEC ( lArnaud 1996 ). The power-law is modulated by the galactic absorption. The 
complete photon spectra for both the observations in the 0.35-10.0 keV energy range could 
not be fitted well with the power-law and galactic absorption model. On the other hand the 
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data and folded model 




Energy (keV) 

Figure 3. Swift/XRT spectrum with best-fit model zwabs*zpowerlaw 



spectra were fitted well by a power-law in the energy range 2.5-10.0 keV. The photon indices 
were found to be 1.35±0.15, 1.37±0.04, 1.38±0.06 and 1.45±0.06 for PN spectrum (Obs-A), 
PN spectrum (Obs-B), MOSl (Obs-B) and M0S2 (Obs-B) spectra respectively. The xV^of 
values are 73.20/66 (PN, Obs-A), 335.49/312 (PN, Obs-B), 274.6/324 (MOSl, Obs-B) and 
388.13/302 (M0S2, Obs-B). The power-law component when extended to lower energies 
(upto 0.35 keV) revealed a soft excess in the 0.35-2.5 keV energy range. The presence of 
soft excess component in the PN and MOS spectra are shown in Fig 1. 
Even though the soft excess in X-rays is seen case of narrow line Seyfert galaxies and it is one 
of the most important components of AGN spectra, its origin is still not well understood. Dif- 
ferent physical models have been used to descri be the soft excess cornponent, such as, fi) th e 



multicolour blackbody radiation fror n the disc (Mitsuda et al. 



ii) thermal Comptoniza tion model (iTitarchuk 



1984 



Makishima et al. 



19861), 



Gie 



rlihski fc Done 



1994J ). (in) smeared wind ab sorption model 



2004J ) and (iv) relativistically blurred reflection model ( iRoss fc Fabian 



20051 ) . Here we use all four models to understand the soft excess component in X-ray spec- 
trum of PMN J0948+0022. The merits and demerits of these models are discussed in the 
next section to choose the best plausible model. 

Before using the physical models, we used phenomenological model broken power-law to 
get an idea of the relative steepness of the soft and hard components of the spectrum. The 
fitting of PN spectrum of Obs-B with a broken power-law yielded photon indices 1.92±0.01 
and 1.46±0.02 below and above the break energy at 1.68±0.02 keV. The x^/dof came out 
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to be 843.7/697. With this estimations of the photon indices we apphed different physical 
models to fit the PN spectrum. 



Multicolour disc blackbody and power-law 

Considering the fact that the soft excess is contributed by the accretion disc of the galaxy 
around the central black hole, we fitted the spectrum with spectral models multicolour disc 
blackbody and a power-law. The complete spectral model used was zwabs*(diskbb+zpow). 
The equivalent hydr ogen column density used to account the galactic absorption was fixed 



at 4.4 X 10^° cm"^ ( Kalberla et al. 



20051 ). The power-law photon index (F) was 1.53±0.02 



indicating a hard power-law and the inner disc temperature was found to be 162±6 eV. The 
spectral fit was reasonably good with x^/dof = 827.45/693 and the fitted spectrum with 
its residual is shown in Figure O The disc temperature obtained here was much higher than 
the inner disc temperature expected for an accretion d isc around a black hole (BH) of mass 



3 X IO^Mq which is the mass of the BH estimated by 



Yuan et al. 



2008h . 



Thermal Comptonization of disc photons 



The other possibility could be that the seed photons originating from the accretion disc 
were Compton boosted by thermal electrons in an optically thick corona to produce the soft 
excess. We have used ZWABS, CompTT - the thermal Comptonization model and a power- 
law component ZPOW to fit the complete spectrum. Considering the fact that the corona 
temperature ~ 0.150 keV (as obtained for inner disc temperature in case of multicolour disc 
blackbody spectrum) and power-law photon index F ~ 1.9 (as obtained for the steeper 



component of the power-law in the phenomenological b roken power-law mode 



optical depth (r) could be estimated from the relation (ISunyaev fc Titarchuk 

1 



) the plasma 



1980h 



TT 



' fcTe 



T + 



3 
2 



(1) 



where kTg is the electron temperature and F is the photon index. This brought out the optical 
depth ~ 52. With the seed photon ter nperature at 30 eV an d the equivalent hydrogen column 
density was fixed at 4.4 x 10^" cm~^ (Kalberla et al.ll2005l ) the complete spectrum was fitted 
using the model zwabs*(comptt-|-zpow). The fitted value of electron temperature was 
O.mto fol keV and that of photon index was 1.55±0.02. The value of x^/dof is 834.67/697. 
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Smeared wind absorption model 



Gierlinski fc Dona (l2004f ) discussed the relativistic smearing of characteristic atomic ab- 
sorption features in a moderately relativistic wind/outflow from the inner accretion disc to 
explain the soft excess in the X-ray spectrum of Seyfert 1 galaxies. The main parameters of 
the model are the absorption column density, ionization parameter (^) and velocity disper- 
sion (cr) of the Gaussian velocity distribution of the wind. Model is described by swiNDl 
within XSPEC. The model swiNDl along with the photo-electric absorption ZWABS and the 
power-law for the background continuum gives jdo] = 809.98/692 with the present data 
from EPIC-pn. The fitted parameters are - the absorption column density : > 25.2 x 10^^ 
cm~^, the ionization parameter : 3.44±0.12 and the velocity dispersion a > 0.5 (in the unit 
of -) and the photon index is 1.71±0.02. 



Relativistically blurred reflection model 



In this model the cold optically thick accretion disc is illuminated by a power-law distributed 
radiation and the accretion disc medium produces a Compton scattered reflected radiation 
along with fluorescence lines. The resultant emission involves the reflected component from 
the disc and the illuminatin g power-law . The relativistic blurring of resultant emission is 



accounted for by convolving a 



Laorl (jl99ll ) profile used to simulate the blurring emission lines 



from an accretion disc around a maximally rotating black hole. The Laor profile estimation 
considers that the accretion disc is geometrically thin, flat and has well defined inner and 



outer radius. The radial depe ndence of disc e r nissiv ity is assume d to be a power- 



model was first introduced by 



Ballantyne et al. 



(12001bf ). 



Ross fc Fabian ( 1993 



Fabian et al. 



) and used bv iBallantvne et al. 



(120021) and 



Ballantyne et al 



updated by including more ionization states a nd more atomic data by 



a w. Thi s 
feoOlal ). 



fl2003h. It was furthe r 



The updated mod e 



and 



Walton et al 



Ross fc Fabian (2005). 



Ross &: FabianI (120051 ) was successfully used by 



(|2012[ ) for a large sample of AGN while by 



Crummv et al. 



Nardini et al. 



fl2006h 



torn for AGN 



Ark 120. The parameters of the model are inner and outer radii of the accretion disc, iron 
abundance of the gas, spectral index of illuminating power-law radiation, emissivity index 
of the disc, ionization parameter and the inclination of the disc. The fitting of the source 
spectrum in the 0.35-10.0 keV energy band requires an additional local absorption due to 
the warm absorber apart from the galactic absorption. The final model used for the fitting 
was wabs*zwabs*kdblur(powerlaw + atable{reflionx.mod}). The equivalent hy- 
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drogen column density was kept fixed at 4.4 x 10^*^ cm~^ (iKalberla et al.l l2005l ) while the 



equivalent hydrogen column density corresponding to the local warm absorber was fitted to 
8.OI2 ^ 10^*^ cm~^. The fitted value of the inner radius of the accretion disc was 1.49]']q^3 Vg 
where Vg is the gravitational radius of the black hole. The inner radius value was very close 
to the inner most stable circular orbit around a maximally rotating Kerr black hole. We 
obtained the lower limit ~ 6 of the emissivity index for the disc. The power component 
came out to harder as compared to the other models used here. The photon index of the 
power-law component was fitted to 1.30 ± 0.06. The reflionx component had disc iron 
abundance 0.211qq5, ionization parameter 2000.0^^947^4^ and illuminating radiation photon 
index 2.70l[J;J^. The fitted spectrum along with the residual is given in Figure |2j The fitted 
parameters for all four models are given in Table [1] 

Even though the above discussions were based on the EPIC-PN spectrum of Obs-B, we also 
fitted the EPIC-PN spectrum of Obs-A and MOSl and M0S2 spectra of Obs-B with the 
above models. The fitted parameter values were generally consistent, but the error bars were 
large. 



Swift/XRT spectrum 

We fitted the XRT data with model zwABS*ZPOW. The galactic absorption and power-law 
model fitted, but due to poor statistics of the data the soft excess in the spectrum could not 
be detected. In fact the fitting was also not of that good quality yielding x^/dof = 62.04/32. 
The fitted spectrum is shown in Figure [3] and the best-fit parameters are given in Table [H 



3.2 Spectral variability 

The lightcurves with 500 s time bin, in the 0.35-2.50 keV and 4.0-10.0 ke V energy bands. 



and th e corresponding hardness ratio are shown in Figure H^. Following IVaughan et al. 



(l2003al ) we also calculated the fractional variability (-F^r) for the lightcurves in the energy 
bands 0.3-0.8 keV at the step of 0.1 keV, 0.8-1.0, 1.0-2.0, 2.0-3.0 and 3.0-5.0 keV. The rms 
spectrum is shown in Figure Hb. The fractional variability did not vary significantly up to 1 
keV and remained around 10%. Beyond 1 keV Fyar showed a trend of increase and reached 
up to 20%. 

We used 100 sec binned lightcurves for the energy ranges 0.3-2.5 keV, 4.0-10.0 keV and for 
the complete energy range 0.3-10.0 keV to calculate the rms amplitude of variability {<Jrms)- 
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Parameters 




Best-fit values 


XMM-Newton (0673730101) : PN 


Model 
r„ (keV) 

^diskbb 


zwabs*(diskbb 4- 
0.162 ±0.006 
44 29+8-09 

1.529 ±0.016 
9.53l[;j? X 10-4 
827.45/693 


ZPOWERLAW) 


Model 
To (keV) 
fcTe (keV) 

^aompTT 
Npl 


zwabs*(compTT 
3.0 X 10-2 

n 1Q7+0.008 

l.lll°;45 X 10-2 
1.55 ± 0.02 

9-93to.29 X 10-* 
834.67/693 


+ ZPOWERLAW) 


Model 

col. dens. (xlO^^ cm^^) 
log 5 

(T (in units of f/c) 


ZWABS*SWINDl*ZPOWERLAW 

> 25.2 
3.44 ±0.12 

> 0.2 

1.71 ± 0.02 

1.77 ±0.04 X 10-3 

809.98/692 



Model 

col. dens. (xl022 cm~ 

Index 

Rin(G) 

Inclination (deg) 

Fe/solar 

Gamma 

Xi 

Npi 
X^l^ 



wabs*zwabs*kdblur(powerlaw+atable(reflionx.mod) 



0.08 



+0.04 



>5.77 



1.49 



+0.21 
-0.13 
,+2.03 



74 22"' 

21+''-^^ 

2 70+''-i'^ 
^•'U_o.o7 

-,+1472 



2000 



947 
+ 1.64 



+0.06 



1 30 

^•'-'"-0.06 
3.21 ± 0.32 X 10-4 

760/692 



Swijt (00031306014) : XRT 


Model 


ZWABS*ZPOWERLAW 




1.58 ±0.11 




2.55 ± 0.17 X 10-3 


X^lv 


62.04/32 



Table 1. The best-fit values of the spectral parameters derived from spectral fitting of X-ray data with different models using 
x2 minimization technique. 



The rms am plitude of variabili ty for a lightcurve with N data points and mean count rate x 
is given by ( lUttley et all 120051 ) 



Or 



1 ^ 



(2) 



The length of each segment of the lightcurve considered was 2000 sec. The variation of 
the calculated rms variability with the mean count rate for the all three energy ranges are 
shown Figure |5l It is evident that, within the limited dynamical range of the count rate, the 
flux vs Grras follows a linear relationship. This is true for the lightcurves in all three energy 
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Bin time: 500.0 s 0.3 




Figure 4. (a) Lightcurve for 0.3-2.5 keV (S), 2.5-10 keV (H) and the corresponding hardness ratio (H/S); (b) rms spectrum 
as obtained from XMM-Newton observation (Obs-B) 

segments. The distribution of points are fitted witli linear equations with slopes 10.53±0.58 
(for 0.3-2.5 keV band) and 10.76±1.41 (for 4-10 keV band). 



4 DISCUSSION 

To understand the spectral behaviour of radio-loud narrow line Seyfert 1 galaxy PMN J0948-I-0022 
in the X-ray energy band we analysed the archival data of the source from XMM-Newton. 
The EPIC-PN spectrum of the source in 0.3-10.0 keV region could not be fitted with a 
simple power-law. The spectrum in the 2.5-10.0 keV energy range fits reasonably well with 
a power-law distribution (x^/dof = 335.5/310) with spectral index 1.37±0.03 while the data 
in the 0.3-2.5 keV energy range show the presence of a distinct soft excess component when 
compared with the power-law extrapolated from the high energy. The soft excess is com- 
monly observed in case of radio-quiet as well as radio-loud narrow line Seyfert 1 galaxies, but 
in the present context it has different importance. Since PMN J0948+0022 is a gamma-ray 
emitter which establishes the existence of powerful jet in the system, so it is not only impor- 
tant to understand the origin of the soft excess in the X-ray domain, it is also very crucial 
to understand its connection with the gamma-ray emission. The latter may shed some light 
on the disc-jet connection which is not well understood in the theory of ac cretion and jet 
formation, even after a lot work have already been done in this area (e.g. 



Neilsen fc Lee 



(120091)). 

Due to the uncertainties in the understanding of the soft excess in the X-ray spectrum of 
NLSl, we considered four different models to fit the EPIC-pn spectrum of PMN J0948+0022. 
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20 I 1 1 1 1 1 1 r 




count rate 

Figure 5. rms-count rate correlation for XMM-Newton data in the frequency range, 5 X 10"* — 5 X 10~^ Hz ; red: 0.3—2.5 keV, 
blue: 4-10 keV POWERLAW component, black: 0.3-10.0 keV. Same coloured lines are the best-fit functions for the respective 
data sets, inset: zoom-in of the POWERLAW component 



Among the models used here it is very difficult to rule out any of them just by considering 
the value of x^/dof. 

The multicolour disc blackbody and power-law model gi ve the inner disc temperature ~ 
0.162±0.006 keV and a power-law index ~ 1.528±0.016. lYuan et all mm estimated the 
black hole mass ~ lO'^-^M© for PMN J0948+0022. For such value of the BH mass the inner 
disc temperature of a standard accretion disc can not be as high as 0.162 keV. This fact 
essentially rules out the possibility that the multicolour disc blackbody can explain the 
observed soft excess. 

In case of thermal Comptonization model, the soft excess in 0.35-2.50 keV energy range 
was modeled by the Comptonization of disc photons in an optically thick cold thermal 
corona. The characteristic photon temperature from the accretion disk was held fixed at 30 
eV which is the typical inner disc temperature for an accretion disc around a black hole of 
mass IO^'^Mq. The fitted value of the corona temperature was found to be O.lSTto oog keV. 
The optical depth of the corona was fixed at 52. The fitting of the power-law component 
yielded a hard power-law with spectral index 1.48±0.03. This hard power-law component 
can be modeled by invoking a hot corona with moderate optical depth (r ~ 1). Such two 



corona scenarios were discus sed in the literature in the context of exp 



spectra of Seyfert 1 g alaxies (IZdziarski et al. 



200 



3 



1996 



Laor et al. 



Magdziarz et al. 



11998 



aining the X-ray 



Dewangan et al 



19971 ). The seed photons for the cold optically thick corona are possibly the 
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photons from the accretion disc while the seed photons for the hot corona can be from the 
accretion disc or photons of the soft excess from the optically thick corona. To distinguish 
the possible origin of the seed photons in the hot corona, we studied the cross-correlation 
between the soft excess and power-law flux. It was found that there was no correlation 
between the two energy bands. If the soft excess photons served as the seed photons in the 
hot corona then there must be a correlation between the bands. Therefore, it was most likely 
that the disc photons of temperature 30 eV served as the seed photons in the Comptonizing 
hot corona and the power-law component was produced in this part of the corona only by the 
process of thermal Comptonization. The main problem of using this model as it was done in 
case of soft excess using the spectral model Comptt was that the observed spectrum with 
XMM-Newton extends up to 10 keV only. It is to noted that the observed spectrum does 
not show any cut-off at this energy. Due to the absence of any cut-off the temperature of the 
corona which determines the cut-off in the photon spectrum could not be constrained. The 
observation of the source in hard X-rays may reveal the thermal cut-off if the power-law was 
indeed produced by the thermal Comptonization. The X-ray flux obtained in the 0.35-10.0 
keV energy range was 4.05 x 10^^^ ergs cm~^ s~^. This led to a luminosity ~ 3 x 10^^ erg 
s~^ which, for a BH mass ~ 1O'''^M0, is O.GLEdd- 

The thermal origin of soft excess leads to a difficulty that the temperatures of the thermal 



electron required to fit the spectrum fall within a very narr ow range ( ~ 0. 



Jin et al 



feoji); 



0.3 keV) ir- 



Done et al 



respec tive of the mass of the accreting black hole. Recently, 
(I2OI2I ) modeled the X-ray spectra of a large number of unobscured type 1 Seyferts using the 
thermal Comptonization model and it was found that the electron temperature varies from 

0.1-0.6 keV while the black hole mass ranges from ~ 10^ IO^Mq. Therefore, the values 

of electron temperature obtained here are generally consistent with same reported in the lit- 



erature. Such a narrow temperature range and its i ndependence on black ho 



that the so ft excess could have had atomic origin (IGierlinski fc Done 



1993 



2OO5I ) . We used the smeared wind absorption model by 



2004 



e mass suggest 



Ross fc 



Gierhnski fc Dond fl2004D to fit 



abian 



the X-ray spectrum of the source. In this model, the power-law component is produced in the 
Comptonizing corona a nd the soft excess is prod uced due to the smeared atomic absorption 
features in soft X-rays ( IGierlinski &: Dondl2004l ). The smearing of absorption features takes 
place in differentially rotating mildly relativistic wind outfiow/jet from the accretion disc. 
The fitting of the spectrum yielded an absorption wind column density > 25 x 10^^ cm~^, 
the ionization parameter is 3.44±0.12 and the dispersion of Gaussian velocity distribution. 
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cr > 0.5. The value of a implied that the outflow was mildly relativistic with Lorentz factor 
^ 1.2. Therefore this model requires the launching of a relativistic outflow to fit the soft 
excess data. 

The relativistically blurred reflection model which also considered the atomic origin of the 
soft excess, gave a reasonably good fit to the observed spectrum. The fitted value the inner 
radius not only implied a rapidly spinning black hole, it also implied that the emission 
region is situated in the innermost portion of the accretion disc. This also implied that 
the reflected radiation would have very blurring due to the gravitational field of the black 
hole. The illuminating radiation had steeper photon index compared to the hard power-law 
component in the spectrum. The hard power-law component might be produced in a corona 
in the outer portion of the disc as mentioned earlier. The estimated luminosity of the source 
came out to be 3.1 x 10'^^ ergs s~^ which is approximately O.GLEdd- 

The other aspect of the present study is the study of rms spectrum. The calculated rms 
spectrum remained almost constant upto 1 keV and it increased beyond 1 keV. The beyond 
1 keV there might be independentl y varying spectral com ponents and the power spectral 
density could be energy dependent (IVaughan et al.ll2003al ). As the error bars are relatively 
large, so to draw any firm conclusion on spectral variability for photon energies beyond 1 
keV longer observations are required. 

The other important aspect of the emission from PMN J0948+0022 is the linear relation 
between the rms variability amplitude and the mean count ra t e. This feature is also seen 



in case of narrow line Seyfert 1 galaxies (Edelson et al. 



from Seyfert 1 galaxies (IVaughan et al. 



2003c 



b|). 



2002 



Uttley et al 



McHardy et al 



2004h and 



(120051 ) showed that the linear 



rms-flux relation is an outcome of a multiplicative process where the amplitude of short 
time-scale variations are modulated by the longer time scale variations. In case of accretion 
flow the long time scale variations are produced at larger radii whereas t he short time scale 
variations are produced at shorter radii of the accretion disc. In fact iLyubarskiil (119971 ) 
proposed a model where the fluctuations originated at the larger radii of the accretion disc 
were considere d to propagate to smaller radii, to explain the observed X-ray variability in 



X-ray binaries. 



Uttley fc McHardyl (120011 ) showed that the linear relationship b etween the 



rms va riability amplitude and flux is a natural outcome of the model proposed by lLyubarskii 

mm . 

Finally, in the 0.1-300 GeV energy range, Fermi/LAT detected 7-ray emission from the 
source. The analysis of the data for one month around the XMM-Newton observation (Obs- 
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B) discussed here, showed a gamma-ray flux of (8.34 ± 0-1 9 ) x -*-PT^ photons cm~^ s" 



Earli er the gamma-ray detection was reported by lAbdo et al.l (j2009bl ): iFoschini et al.l (12010 



20121 ). The multiwavelength spectral data were modeled by using synchrotron process for the 



radio-to-optical band; X-ray and 7-ray spectra were fitted with self-synchrotron Compton 
and external Compton processes respectively. The radiation is produced in a region in rel- 
ativistically moving jet. In those works the X-ray spectra were considered from Swift/XKT 
in the 0.2-10 keV band and the XRT spectra were fitted with simple power-law. But in the 
present work it is shown that the X-ray spectrum in the 0.2-10 keV energy band was not 
a power-law. The near simultaneous Swift/XBT data were analysed here and it was shown 
that the data could not reveal the soft excess component due to poor statistics. Therefore 
the consideration that the X-ray spectrum was a simple power-law and originated in a rel- 
ativistic jet might be misleading. It indeed has a soft excess component from the accretion 
disc. This must be taken into account while modelling the multiwavelength spectral energy 
distribution. It is also important to understand the relationship between the soft excess 
component and the very high energy component of the spectrum which is beyond the scope 
of this paper. 



5 CONCLUSIONS 

We have analysed the archival X-ray data from XMM-Newton , Swift/XBT and 7-ray data 
from Fermi/LAT for the radio loud narrow line Seyfert 1 galaxy PMN J0948-I-0022. This 
conclude the following. 

• The spectral analysis of X-ray spectrum from XMM-Newton establishes the presence 
of soft excess in the 0.3-2.5 keV energy band and a hard power-law beyond 2.5 keV. 

• Due to the uncertainties in the understanding of the soft excess in the X-ray spectrum 
of NLSl the complete spectrum is fitted with four different spectral models - (1) multicolour 
black body and power-law, (2) thermal Comptonization and power-law, (3) smeared absorp- 
tion wind model and (4) the relativistically blurred reflection model. The models (1) and (3) 
are completely ruled out, but models (2) and (4) are statistically comparable to each other. 
Due to the fact that the temperature of the optically thick corona where the soft excess is 
produced by thermal Comptonization, falls in a very narrow range irrespective of the mass 
of the central black hole, the relativistically blurred reflection model is preferred. 

• The simultaneous multi-wavelength fitting of the spectrum should take care of the soft 
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excess part of the spectrum properly. The power-law component in the X-ray spectrum 
might also be originated in the accretion disc corona. The SSC process might not be the 
correct description of X-ray emission in 0.3-10.0 keV energy range. It is also important to 
understand the possible connection between the soft excess and the gamma-ray emission in 
the source. 

• It is important to study the variabihty properties of this source in detail. This study 
will be helpful to understand the underlying radiation emission process and its impact on 
the variability of the source. 
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